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I. INTRODUCTION 

The sustained interest in Gross-Neveu (GN) models in 
1-1-1 dimensions [ij stems to a large extent from their chi- 
ral properties. Thus the simplest model with Lagrangian 

c = v^ii'^a^^ + \g\4>^f (1) 

(suppressing flavor indices, i.e., "ifjip = X^^i V'fcV'fe stc.) 
has a discrete chiral Z2 symmetry 

i' -> 15^^, (2) 

whereas the chiral GN model or, equivalently, the two- 
dimensional Nambu-Jona-Lasinio model (NJL2) 0, 

£ = ^ird^i; + + \g\i^\i,i^)\ (3) 

possesses a continuous chiral U(l) symmetry, 

?/> ^ e'"''^?/>. (4) 

Chiral symmetry and in particular its breakdown man- 
ifest themselves in such diverse physical phenomena as 
dynamical fermion masses, the meson spectrum, topo- 
logical effects in the structure of baryons, and rich phase 
diagrams at finite density and temperature with various 
types of homogeneous and solitonic crystal phases, see 
the introductory review article 0] as well as the recent 
updates in By adding a bare mass term to the 

Lagrangian, one breaks the chiral symmetry explicitly 
and gets additional insights into the symmetry aspects 
of both models 0, Q. Nevertheless, studies of models 
(dl) and ([3]) with their strikingly different properties have 
remained somewhat disconnected. 

In the present work, we propose and solve a simple 
field theoretical model which interpolates continuously 
between the Lagrangians H]) and ([3]). Our motivation 
is to get a better understanding of how the conspicuous 
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differences in the phase diagrams and baryon structure 
come about. Moreover, we would like to explore an al- 
ternative mechanism for breaking chiral symmetry ex- 
plicitly, different from the usual bare mass term. To this 
end, we consider a Lagrangian similar to Eq. ^ , but with 
different (attractive) scalar and pseudoscalar couplings, 

C = {pird^yj + ^g^i'^f + ^G^^ij5^f. (5) 

By varying from to g^, we generate a family of theo- 
ries interpolating between the GN and the NJL2 models. 
The idea to generalize the GN model in this fashion is not 
new. Thus for instance, Klimenko has studied a closely 
related problem long time ago However, since the 

role of inhomogcneous condensates has only been appre- 
ciated in recent years, there is almost no overlap between 
the present work and these earlier studies. 

The methods which we shall use in our investigation 
have been developed during the last few years in an effort 
to clarify the phase structure of massless and massive GN 
models. As a result, we have now at our disposal a whole 
toolbox of analytical and numerical instruments. The 
most important keywords are: the derivative expansion, 
asymptotic expansions, perturbation theory, Ginzburg- 
Landau (GL) theory and numerical Hartree-Fock (HF) 
approach including the Dirac sea. This will enable us to 
solve the generalized GN model ([5]) in a rather straight- 
forward fashion, although the model is far from trivial. 
Its two limiting cases, the standard massless GN and 
NJL2 models, can both be solved analytically. This is 
unfortunately not true for the generalized model which 
in this respect is closer to the massive NJL2 model Q. 

This paper is organized as follows. We present our 
computations and results starting with mostly analyti- 
cal work and ending with purely numerical results. The 
logic of the HF approach demands that we begin with a 
discussion of the vacuum, dynamical fermion mass and 
coupling constant rcnormalization in Sec. |TT1 Sec. IIIII is 
dedicated to fermion-fermion bound states (mesons) and 
scattering. In Sec. IIVI we solve the theory in the baryon 
sector as well as for low density soliton crystals in the 
vicinity of the chiral limit, using a kind of chiral perturba- 
tion theory obtained from the derivative expansion. We 
then begin our study of thermodynamics at finite tem- 
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peraturc and chemical potential with an investigation of 
the tricritical behavior near the chiral limit in Sec. [V] In 
Sec. IVII the microscopic GL approach underlying Sec. |V] 
is extended to more general coupling constants, and the 
tricritical point of the generalized GN model is deter- 
mined exactly. Some technical details are deferred to 
the appendix. Sec. I VIII is devoted to the full phase dia- 
gram of the generalized GN model for arbitrary coupling 
constants, chemical potential and temperature, only ac- 
cessible via a numerical relativistic HF calculation. As a 
by-product, we also present information about baryons 
away from the chiral limit. The paper ends with a con- 
cluding section. Sec. I Villi 



II. VACUUM, DYNAMICAL FERMION MASS, 
RENORMALIZATION 



Consider the Lagrangian of the generalized GN model 
with two coupling constants in 1-1-1 dimensions, Eq. ([5|). 
For = g^, it coincides with the one from the mass- 
less NJL2 model, Eq. (jSj. For = 0, we recover the 
massless GN model, Eq. ([T]). The case G^ > can be 
mapped onto G^ < g^ by means of a chiral rotation about 
a quarter of a circle. 



^i75 7r/4 



(6) 



Since this is a canonical transformation, we may assume 
< G^ < g^ without loss of generality. Hence the gener- 
alized GN model can serve as a continuous interpolation 
between two well-studied model field theories with dis- 
tinct symmetry properties. Notice that the generalized 
Lagrangian ([S|) always has the discrete chiral symmetry 
■0 ~^ 75^ under which iptp and V'^TsV' change sign. The 
continuous chiral symmetry ip — > e'"'''^ ip is only recovered 
at the point g"^ = G^. 

To find the vacuum in the large N limit, we introduce 
homogeneous scalar and pseudoscalar condensates, 



M = -02(^175^-). 
The Dirac-Hartree-Fock equation 

(-75i5x + -f°m + h'^M) ^^Eil; 
then yields the single particle energies 

E = zbv/fc^ + m2 + Af2 
and the (cutoff regularized) vacuum energy, 

p — _ 



(7) 
(8) 
(9) 



dfc 



A/2 



27r 



A^ m? -t- 

-8^ + 

,2 



47r 

'2Ng^ 27VG2' 



In 



2Ng^ 



2NG^ 



A2 



- 1 



TO 



(10) 



If we choose the following relations between the UV cutoff 
A/2 and the bare coupling constants g^,G^, 



Ng^ 



In A 
In A 



6, 
6, 



(11) 



^vaci'm-, M) is well defined in the limit A — > 00 (dropping 
the irrelevant quadratic divergence) and given by 



M 



[ln(TO^-,M^)-l]+^ + ^. 

(12) 



Minimize fvac with respect to to, M , 
TO [2^1 + In (n 



A/ [2^2+ In I 



(13) 



These equations only admit a solution with nonvanishing 
TO and M if ^1 = ^2 = -^\n{m'^ + M^). This takes us 
back to the NJL2 model with its infinitely degenerate 
vacua along the chiral circle of radius + IvP . The 
other options are to 7^ 0, M = 0, ^2 unspecified and 



1 to 
or else m = 0, M 7^ 0, ^1 unspecified and 



(14) 



(15) 



The vacuum energy is lowest for m 7^ if ^1 < ^2 and for 
M 7^ if ^1 > ^2 • In view of the remark below Eq. ([6|) , we 
may adopt the first scenario. Choosing units such that 
TO = 1 and denoting ^2(> 0) by ^ from now on, we finally 
get the renormalization conditions (gap equations) 



Ng^ 

TT 

7VG2 



= In A, 



Ng^ 



(16) 



With the help of these relations, all physical quantities 
can be expressed in terms of the scale to (set equal to 
1) and the dimensionless parameter ^ which serves to 
interpolate between the massless NJL2 {£, = 0) and GN 
(^ — 00) models. This expectation is borne out in the fol- 
lowing sections, supporting our renormalization method. 



III. MESON SPECTRUM AND 
FERMION-ANTIFERMION SCATTERING 

In the large N limit, fermion-antifermion bound and 
scattering states can conveniently be derived via the rel- 
ativistic random phase approximation (RPA) [ll|, [l^ . 
Since the scalar and pseudoscalar channels decouple and 
the HF vacuum is the same as in the GN or NJL2 model, 
this analysis requires only minor changes of the standard 
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calculation for the NJL2 model. Consider first the bound 
state problem. The scalar channel has been spelled out in 
all detail in Ref. where it is shown that the eigenvalue 
equation assumes the form 



1 = 2Ng' 



dk 



u{k)v{k - P)u{k - P)v{k) 



2tt 

E{k - P, k) 



'£'^{P)-E^{k~P,k) 



(17) 



Here, P is the total momentum of the fermion- 
antifermion system, w, v are positive and negative energy 
HF spinors, and 



E{k) = Vfc^ + 1, E{k',k)^E{k')+E{k). (18) 

The energy of the meson is denoted by £{P) = 
+ A^^. An analogous computation in the pseu- 
doscalar channel gives 



1 



dk 

-2NG'^ I —u{k)i-f5v{k ~ P)u{k ~ P)ij5v{k) 
2tt 



E{k - P, k) 



£^{P)-E^{k-P,k)' 
Use of the identities 

u(k)v{k - P)u{k - P)v(k) 



(19) 



4 + p2 -E^{k- P,k) 



AE{k)E{k - P) 



(20) 



~L(k)'v^Qv{k — P)u(k — P)i75t;(fc) = — 



P"^ - E^{k- P,k) 



AE{k)E{k~ P) 
(21) 

puts these eigenvalue equations into the more convenient 
form 



1 



1 = 



dk 



1 



1 



2 J 271 \E{k - P) E{k) 
4 + P^ - E^{k - P,k) 
^ £^{P)-E^{k-P,k) ' 
NG^ f dk f 1 1 



2 J 271 \E{k - P) E{k) 

P^ - E^{k- P,k) 
£^{P)-E^{k-P,k)' 



(22) 



If we regularize the momentum integrals with the same 
cutoff A/2 as used in the treatment of the vacuum energy 
and use the renormalization conditions Eqs. (jl6p . we get 
the renornialized eigenvalue conditions 



= 



TT 



dfc 



27r \E{k - P) E{k) 

£^{P)~E^{k~P,k)' 
dk / 1 1 



27r \E{k - P) E{k) 
p2-£2(p) 

'£'^{P)-E^{k-P,k)' 



(23) 



(24) 



now free of divergences. Eq. (j23p is the same as in the 
standard GN and NJL2 models and gives the familiar 
result for the scalar (a) meson mass, M = 2. The right- 
hand side of Eq. is independent of P and can readily 
be evaluated in the cm frame of the meson (P = 0), 



2 
1 



e = / dfc 



1 



Vfc2+T(7W2_4_4fc2) 



: arctan ■ 



1 



with 



4 



(25) 



(26) 



Solving the transcendental equation (j25p numerically, the 
pseudoscalar (tt) meson mass is found to rise from = 
at ^ = to 2 at ^ — >■ 00, see Fig. [T] The first limit is as 
expected - this is the would-be Goldstone boson of the 
NJL2 model. The 2nd one is surprising at first glance, 
since we are supposed to reach the GN model in this 
limit. The GN model does not have any pseudoscalar 
fermion-antifermion interaction, let alone a bound state. 

To better understand what is going on, we briefly turn 
to the fermion-antifermion scattering problem. Since the 
RPA equations have a separable kernel with one-term 
separable potentials in the scalar and pseudoscalar chan- 
nels, this is straightforward [isj . The energy dependence 
of the scattering matrix is encoded in the following func- 
tions of the Mandelstam variable s. 



1 ^ A/«2 f dfc 1 4fc^ 

+ J 27r vT+F" s-4(l + fe^)+ie 

1 I ATr'2 r dfc 1 4(l + fc2) 

^ + ^^^ J 2¥7T^s-4(l+fc-^)+ie 



(27) 



Upon isolating the divergent part of the integrals and 
using the renormalization conditions, this becomes 



7T ZTT 



lis) 



dfc- 



1 



1 



(28) 



Vl + fc2 s - 4(1 + fc-2) + ie 
where the integral I{s) can be evaluated in closed form. 



lis) = - 



lis) 



In- 



: arctan 



4- s 



is < 4) (29) 



(s>4) 



(30) 

To- has the expected pole at s = 4 corresponding to the 
marginally bound scalar meson with A4 = 2. The pole 
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FIG. 1: Masses of a and tt mesons vs. ^ in the large A'^ limit 
of the generalized GN model, obtained from Eqs. (|23j|26p . 



of in turn coincides with the mass of the pseudoscalar 
meson, see Eqs. (|25I26|) . According to the 2nd hne of 
Eq. p8)) . the strength of the pseudoscalar scattering ma- 
trix vanishes hke ~ 1/^ for ^ — > cxd. We therefore arrive 
at the following picture: As ^ — >■ oo, the pseudoscalar in- 
teraction vanishes, in accordance with the expected GN 
limit. However, since an arbitrary weak attractive in- 
teraction is sufficient to support a bound state in 1-1-1 
dimensions, the pseudoscalar bound state pole persists, 
the binding energy going to zero. As we shall see later 
on, this decoupled tt meson has no influence on any other 
observables of the model in the large N limit, so that it 
does not really upset our goal of interpolating between 
the NJL2 and GN models. 



IV. BARYONS AND SOLITON CRYSTALS AT 
SMALL C AND LOW DENSITY 

The derivative expansion is a standard technique 
to deal with quantum mechanical particles subject to 
smooth potentials [l3, EH- In Ref. [l^ it has been 
adapted to the particular needs of the HF approach 
for low dimensional fermion field theories. In effect, it 
amounts to integrating out the fermions in favor of an 
effective bosonic field theory, where the scalar and pseu- 
doscalar fields can be identified with the HF potentials re- 
lated to the composite fermion operators ipip and ipij^ip. 
For baryons in the massive NJL2 model it leads to a chiral 
expansion in closed analytical form [l^. Note that this 
method can only handle fully occupied valence levels at 
present. 

Since the HF equation in the problem at hand has the 
same form as in the NJL2 model, we can take over the 
derivation of the effective action from Ref. [l^ almost lit- 
erally. The Dirac-HF equation is written as in Eqs. (|7l8p 
except that the scalar (S) and pseudoscalar (P) conden- 



sates in the baryon state are x dependent, 

[-I5'^d, + -/°S{x) + i-f^Pix)] i> = Eij, (31) 

with 



S = -g'{H'), 
P = -G2(^i75V>. 



(32) 



As is well known, the HF energy can be written as the 
sum over single particle energies of occupied orbits and 
a double counting correction. Only this last part is dif- 
ferent in the present case. Due to the renormalization 
condition (|16p . it depends on the parameter 5, 
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p2 



s- 



2Ng^ 2NG^ 



I p2 c 

_1 inA+^P^ (33) 

2n 2tt ^ ' 



The cutoff dependent term cancels exactly the logarith- 
mic divergence in the sum over single particle energies. 
Only the last term in Eq. ([33| is different from what it 
was before. Consequently, we can simply take over the 
effective action from Ref. [l^, set the confinement pa- 
rameter 7 = (vanishing bare fermion mass) and add 
the new contribution proportional to ^ from Eq. ()33|) . 
Adopting polar coordinates in field space. 



5 - iP = (1 + A)e 



2ix 



(34) 



and working at the same order in the derivative expansion 
as in [1^ , we then get at once the energy density (' = 
and denotes the 4th derivative of x) 



27rf 



= e(l + A)2 sin2(2x) + - \U'f + ^(x'")' 



— (y 
140^^ 



Y - —ir 



+ A^ + -(A')^ 



— {X"Y - -XiX'f - — A* + -Xixy 
120^ ' 6 ^ ' 12 3 ' 



5 2 



(35) 



We have to vary the energy functional with respect to 
A and x solve the Euler-Lagrange equations, then 
compute baryon number and baryon mass. Although we 
shall follow the same procedure as in Ref. [l^ , the results 
will be quite different, reflecting the different ways in 
which chiral symmetry is broken in these two models. For 
simplicity, take first the case of the leading order (LO) 
derivative expansion. Here, we only keep two terms in 
the energy density. 



2^f =Csin2(2x) + (x')' 



(36) 



Rescaling the chiral phase field and its spatial argument 
as follows. 



1 



X{x) = -e{y), y = 2^2 



(37) 



we recognize the (static) sine-Gordon action (' = dy) 



An 



-£ 



COS0 + 1. 



(38) 
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The Eulcr-Lagrangc equation is the time-independent 
sine- Gordon equation 



= smf 



(39) 



so that the baryon can be identified with the sine-Gordon 
kink 



9 — 4 arctane'' 



(40) 



But unhke in the massive NJL2 model, this object has 
baryon number 1/2, exactly like the kink in the standard 
GN model (with fully occupied zero-mode), 



dx 



X 



1 



[x(oo) - x(-oo) 



(41) 



Here we have used the topological relationship between 
baryon number and winding number of the chiral phase 
pH [l6j . The mass of this kink-like baryon is found to be 



Mb 



TT 



rrhr_ 
2tt 



(42) 



where, in the 2nd step, we have made use of Eq. ([25|) to 
LO in ^ and denoted the pion mass by rriT^. 

In the same vein, higher order calculations closely fol- 
low Ref. [1^. We find it useful to switch from the pa- 
rameter ^ to ttItt by means of Eq. (j25p , 



1 



1 



1 



and to expand x ^-nd A into Taylor series in mjr. 



X 
A 



Xo 



-'71^X2 



'71^X3, 



m^Ai + TO^A2 + m^Aa. 



(43) 



(44) 



The Euler-Lagrange equations corresponding to the ef- 
fective action ([55]) can then be solved analytically with 
the NNNLO resuhs {y = m^x) 



Xo 



arctan e 
1 



y 



Xi 



4 cosh^ y 
1 sinh y 



16 cosh^ 



V 



A2 = - 



1 lOcosh^y- 13 



96 



cosh^ y 



(45) 



X2 



1 sinh y(ll cosli^ y — 26) 



2304 



cosh'* y 



X3 



1 562 cosh" y - 3090 cosr y -I- 281 1 
^5760 cosh^ y 

sinh y (6271 cosh'' y + 29588 cosh^ y - 26784) 



1382400 
The baryon mass becomes 



cosh^ y 



Mb 
N 



2tt 



1 ? 



13 
3600' 



1193 
705600 ' 



(46) 



As the whole winding number of x resides in the LO term 
XOj baryon number is always 1/2. Therefore the complex 
potential S ^ iP traces out half a turn around the chiral 
circle. This is confirmed by plotting S and P, showing 
kink-like behavior of S like in the massless GN model, 
see Fig. 2. The presence of a non- vanishing P signals 
that we are dealing with a new kind of solitonic baryon 
here which did not show up yet in any other variant of 
the GN model family. 

Let us now turn to periodic solutions of the Euler- 
Lagrange equations in the derivative expansion. They 
are expected to approximate systematically the ground 
state of matter at low densities and in the vicinity of the 
chiral limit ^ = 0. Since the resulting expressions are 
rather lengthy, we only give them up to NNLO here. 



Xo 
Ai 

Xi 

A2 

X2 



TT 1 

— I — am 
4 2 

1 2 

4 



1 + 1 
24 16 

13 C 
96 " 24 

c 



c 



sn 



24k2 

1 1 
24 



dnZ 



24k2 



sn 



96^2 



24k2 



sn cn dn Z 

(«2 „ 



V 576^2 576k2 
59k2 - 44' 



2304^2 

V 576^4 288^4 

13 C_ 
1152 ^ 96 ^ 576 



(6H-30k2)C 
2880^2 



(61 4-30k2)C 
2880^4 

'2 



dnZ 



C r,2 

^sn cnZ 

576^2 



288^2 96^2 



dn sn^ Z 



Here, 



K 

E' 



(47) 



(48) 



E, K are complete elliptic integrals of k and am, sn, cn, dn 
and Z are standard Jacobi elliptic functions with spatial 
argument 

z = —X (49) 

and elliptic modulus k. The mean density can be simply 
inferred from the period of the crystal. 



P = 



4kK 



(50) 



By way of example, we show in Fig. [3] the scalar and 
pseudoscalar potentials corresponding to ^ = 0.2 (as in 
Fig. [2]) and the density p = 0.05. Again the convergence 
seems to be very good. 
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FIG. 2: Scalar {S) and pseudoscalar (P) potentials for 
baryon in the derivative expansion, ^ = 0.2, m-n « 0.8389. 
Dashed curves: LO (sine-Gordon), solid curves: NNNLO, see 
Eqs. g5]). 
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FIG. 3: Soliton crystal for generalized GN model, ^ — 
0.2, rriTT ~ 0.8389, p = 0.05. Dashed curves: LO (sine Gor- 
don), solid curves: NNLO, see Eqs. (fi?!) . 



us with a way of testing the results from the derivative 
expansion. Starting point is the divergence of the axial 
current in the generalized GN model 

« = -2(.g2-G2)^V^i75V' 



2 {S^pi-y^^ - Pip^p) 



2N^ 



SP, 



SP 



(52) 



where we have taken a ground state expectation value 
and used large factorization. Owing to the properties 



7° 



J,=f (53) 
specific for 1-1-1 dimensions, we get for stationary states 

2A^^ 



dip{x) 



-S{x)P{x). 



(54) 



Twofold integration for the baryon case then leads to a 
sum rule relating baryon number directly to an integral 
over the HF potentials S, P, 



2N£, 



2C 



dx'S{x')P{x') (55) 
dx I dx'S{x')P{x') 

(56) 



^ / dxxS{x)P{x) 



In the last step, partial integration was used. Inserting 
the results for S,P from the baryon, i.e., 



S = +(l + A)cos(2x), 
P = -(1 + A)sin(2x), 



(57) 



with X, A from Eqs. (P5|) . we find that the sum rule is 
only violated at 0(m|). This is a good independent test 
of a considerable amount of algebra behind the derivative 
expansion. 



Since the derivative expansion is anyway expected to 
be most useful at low densities, we note the following 
simplification in the low density limit: for k — > 1, we 
can use the approximation ^ w and keep k only in the 
arguments of the Jacobi elliptic functions. Expressions 
(j47p then reduce to periodic extensions of the baryon 
results obtained by simply replacing 

1 sn(z,K) 

coshy — r, smhy — > — -, r (51) 

cn(z,K) cn(z,K) 

in Eqs. P5)) . 

Finally, we derive a sum rule for the baryon number 
of a single baryon, following Ref. [l7| . This will equip 



V. PHASE DIAGRAM NEAR THE NJLz 
TRICRITICAL POINT = 0) 

We start our investigation of the phase diagram of the 
generalized GN model by zooming in onto the tricritical 
point at ^ = 0, i.e., of the NJL2 model. In Ref. [l^ it 
was shown that this region is well suited for the derivative 
expansion, which here leads to a (microscopic) Ginzburg- 
Landau type theory. In that work, chiral symmetry was 
broken as usual by means of a bare fermion mass term. 
Here instead we break it by choosing two slightly dif- 
ferent coupling constants in the scalar and pseudoscalar 
channels. The central quantity of interest is the grand 
canonical potential which differs in these two cases only 
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by the double counting correction. Since the latter is 
independent of temperature and chemical potential, the 
situation is very similar to the one in the preceding sec- 
tion. Once again we can take over the effective action 
from the literature about the massive NJL2 model p^ . 
The only necessary modification is to replace the double 
counting correction term coming from the bare mass by 
the one proportional to ^, cf. Eq. For the present 

purpose, it is advantageous to combine the HF potentials 
S, P into one complex field = 5— '\P. The result for 
the grand canonical potential density to the order needed 
here (dropping a field independent part) then becomes 



with 



a2|0| 



2 I rx. 



a2 = 

= 

04 — 



/*) + a4(|</'|^ + |0?) + |^(5</ 



2tt 



[ln(47rr) + 3?^'(z)] 



87r2r ^ ^ 



647r3r2 



and 



f 

z = — 



(59) 



(60) 



2 27rr 

We denote the digamma and polygamma functions as 

d d" 
«'(z) = — Inr(z), ^-("'(2) = — *(z). (61) 

In the chiral limit (^ = 0), the tricritical point is located 
at 



with the constant 

16e2C 

a = ~ 6.03198, (66) 

the reduced grand canonical potential density 

^-cff = ^^-cff (67) 
becomes indeed independent of ^, 



(68) 



(58) The Euler-Lagrange equation 



- 2ivip + (cr^ - v^)ip - 2|(^|V (69) 

differs from the complex non-linear Schrodinger equation 
by the term ^ (p* . This has prevented us from finding 
the solution in closed analytical form. Let us first de- 
termine the expected 2nd order phase boundaries. The 
phase boundary between massless and massive homoge- 
neous phases can easily be found by minimizing 5" off with 
the ansatz ip = m and setting m = in the condition for 
the non-trivial solution. The result in the new coordi- 
nates is the straight line 



a — i^. 



(70) 



Next consider the phase boundary separating the crystal 
phase from the chirally restored (m = 0) homogeneous 
phase. Here we use the ansatz (see Sec. IV of Ref. [l^ 
for the justification) 



= 0, Tt 



Tc = — 

TT 



(62) 



with Euler's constant C « 0.577216. Following Ref. [18|, 
we expand the coefficients ([59| of the GL effective action 



around the tricritical point ()62p. 




^ f C(3)e-^^A.^ 
on 


as s 









(63) 



with r = \/Tc — T. The ^-dependence can now be re- 
moved as follows. Rescaling the field and the coordiante 
according to 

(j){x) = ^^/^V3(u), u^^^^^x 
0'(x) = e^(u), 0"(x) = e^^0{u) (64) 

and introducing rescaled thermodynamic variables 

2/x fTT T 



^1/2 



Te ^1/2 



(65) 



ip = Co cos{qu) + ido sin((?u) 



(71) 



and evaluate the spatial average of ^off, keeping only 
terms up to 2nd order in cq, doi 



(#cff) = Miicl + 2Mv2Codo + M22dl, (72) 



with 



M12 = -vq 

M22 = ]-{a + q^ + -cj'^). 



(73) 



As explained in Ref. |8|, the phase boundary is now de- 
fined by the conditions 



d 

detM = 0, — detA^ = 0, 
oq'^ 



yielding the critical curve 



^a(8i/2 - a) 
4l^ ' 



(74) 



(75) 
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The wave number q obeys 



(76) 



The tricritical point can be identified with the point of 
intersection of the two critical curves (|70p and ([75]). 



at = i^t 



(77) 



Going back to the original, unsealed variables, this trans- 
lates into 

1 



(78) 



Notice that q vanishes at the tricritical point. We ex- 
pect that a third critical line ends at the tricritical point, 
namely the 1st order phase boundary separating the crys- 
tal from the massive Fermi gas phase. It has to be de- 
termined numerically. To this end, we insert the Fourier 
series ansatz 



= ^ c„ cos[(27i + l)qu] +i^dn sin[(2n + l)qu] (79) 



into Eq. (j68p and minimize the effective action with re- 
spect to the parameters Cn,dn and q. By keeping only 
wave numbers which are odd multiples of q, we restrict 
ourselves to potentials which are antiperiodic over half a 
period, 



ip{u + n/q) = -ip{u). 



(80) 



This kind of shape is indeed favored by the minimization, 
as was the case for the massless GN model. It shows that 
discrete chiral symmetry and translational symmetry are 
broken down to a discrete combination of the 2 transfor- 
mations, namely 



(81) 



from which Eq. (|80p for bilincars follows. In practice, 
we found that it is sufficient to keep co,ci,do,di in the 
expansion ([7^ . Comparing the reduced grand potential 
with the one from the homogeneous massive solution, we 
can locate the phase boundary. The result of the calcu- 
lation is shown in Fig. Altogether with the two 2nd order 
phase boundaries discussed above. Due to the rescalings, 
this is a kind of universal phase diagram which contains 
all information about the actual phase diagram in the 
vicinity of the tricritical point at ^ = 0. By undoing the 
rescaling we can reconstruct the phase diagrams for small 
^ values in a limited region of the (/i, T) plane. This is 
shown in Fig. 5. Here one sees nicely the transition from 
the behavior qualitatively familiar from the GN model 
to the one from the massless NJL2 model. The angle 
between the two phase boundaries delimiting the crystal 
at the tricritical point is consistent with zero, just like in 
the standard GN model. 



a 




crystal 



FIG. 4: Rescaled phase diagram near the tricritical point of 
the NJL2 model. Straight line: 2nd order phase boundary, 
Eq. (|70p . Dashed curve: 2nd order phase boundary, Eq. (|75|l . 
Solid curve: 1st order phase boundary, numerical calculation. 
The 3 critical curves meet at the tricritical point at = ft = 
y^/2. The parameter ^ has been eliminated by the choice of 
variables, see Eq. (fSSj) . 




FIG. 5: Reconstructed phase diagram of generalized GN 
model near the tricritical point of NJL2 model for ^ — 0.0001, 
0.0002, 0.0004, 0.0007, 0.001, 0.002, 0.004, 0.007, 0.01, from 
left to right. All curves are obtained from the ones shown in 
Fig. m but ly, T values up to w 50 are needed for the smallest 
^ value. 



VI. EXACT TRICRITICAL BEHAVIOR FROM 
GINZBURG-LANDAU THEORY 

As ^ varies from to 00, the tricritical point of the 
generalized GN model moves from the NJL2 to the GN 
tricritical point, i.e. from fi = 0,T = 0.5669 to fi = 
0.6082, T = 0.3183. Since the HE potential 4, = S - iP 
vanishes at the tricritical point and its period is expected 
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to diverge, the derivative expansion should be sufBcicnt 
to determine the exact tricritical behavior for all ^. As a 
matter of fact, this will enable us to determine analyti- 
cally the location of the tricritical point as a function of 

We will also be interested in the behavior of the phase 
boundaries in the vicinity of the tricritical point. It turns 
out that the region of validity of the GL theory as defined 
in Eq. (j58p shrinks rapidly with increasing ^. One of the 
reasons is the fact that both a2 and vanish at the GN 
tricritical point, so that it would be necessary to go to 
higher orders in the derivative expansion for large ^. To 
keep the analytical work reasonably simple, we therefore 
analyze the phase boundaries only for moderate ^ values. 

We start once again from the GL effective action . 
Consider first the homogeneous phases. The constant 
ansatz (j) ~ m yields 



cff 



a2Tn 



aim. 



(82) 



Minimizing with respect to m, we find either m = or 



(«2 < 0). 



(83) 



We thus recover the well known result for the phase 
boundary between massless and massive Fermi gas 
phases, namely 



or, parametrically (parameter v), 



(84) 



47r 



z = - +\- 



2ti 



DT. 



(85) 



Next consider the 2nd order phase boundary between 
crystal and massless homogeneous phase. As in Sec. |Vl 
the ansatz 



or, equivalently. 



27ra4 




= 



Q!2 

a4 



2al 



(90) 



These two equations determine Q and the critical curve 
in the (/^, T) plane. The tricritical point must lie on this 
curve and on the curve a2 = 0. This gives the conditions 
Q = and 



27ra| 



a4 



(91) 



where the right hand side is to be evaluated at the tri- 
critical point. Using Eqs. ([5^ . we finally arrive at the 
following parametric representation of the dependence of 
the tricritical point {fit,Tt) on ^ (parameter i>t), 



2 [3^1/(1) (z,)]' 



/.tt = i>tTt. 



3fi*(2)(zt) 



Zt 



2tt 



(92) 



This result should hold exactly in the generalized GN 
model, since GL theory becomes rigorous at the tricritical 
point. It has the correct limits for ^ — )• (NJL2) and 
^ — > 00 (GN), as follows immediately from the vanishing 
of as and 04, respectively. Moreover, by expanding in vt 
we recover the asymptotic behavior of {fj.t,Tt) for ^ — > 
found in Sec. El cf. Eq. 

We now determine the shape of the phase boundaries 
near the tricritical point for finite ^ values. To this end, 
we measure chemical potential and temperature from the 
tricritical point (at fixed ^), 



Co cos{Qx) + ido sm{Qx) 



(86) 



is adequate for a continuous phase transition which can 
be treated in perturbation theory. The spatial average 
of the effective action, keeping only quadratic terms in 
(co,do)i then becomes 



where 



cff/ = 

Mn 
M22 



Miicl -f 2>li2Codo + M22dl 



(87) 



{a2 + aiQ 



1 



- I a2 + oiiQ 



27r 



The 2nd order phase boundary is again defined by 

d 



ActM = 0, 



dQ- 



■ dctM = 



(88) 



(89) 



T 



Tt 



(93) 



We then rotate the coordinate frame in the {S, r) plane 
such that the new axes are tangential and normal to the 
homogeneous phase boundary 02 = 0, 



cos f — sm f 
sin 9 cos t 



with 



smf = — 



cos ^ = ^ 



(94) 



(95) 



We have defined 



C = 27r + i>t/i, n^^ll + C^, h^^^('>izt). (96) 

Due to the cusp, the phase boundaries lie in the region 
around the tricritical point where 



(97) 
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In this region, the Taylor expansion 

012 = a22£^ + ■•■ 

"3 = 030 + a3i£ + 0326^ 
2 

a4 = + a4i£ + a42e 



The choice 



(98) 



holds with calculable coefficients given in the appendix. 
We first determine the shape of the 2nd order phase 
boundary from (|90l9ip and To leading order in 

e, we find 



030(2031040 — 030041) 



2olo 



(99) 



and the following condition for the phase boundary, 

= 4a3oa3ia4oa4i + 4a22a|o - -^al^al^ ~ agoal^. (100) 

We can also determine the ratio (Iq/cq of imaginary to 
real amplitudes. 



do 

Co 



2031040 — 030041 



2030O40 



(101) 



Computing the 1st order phase boundary is the most 
complicated task. Let us decompose (j) into real and imag- 
inary parts and assume the following LO behavior in e, 

4> = F + iG 

y = e'/^x 

Fix) = eFo(y), F'{x) ^ e^/^Foiy) 

G{x) ^ e^/^Goiy), G'{x) = e^Go{y) (102) 

These assumptions will be justified a posteriori once we 
have constructed a consistent solution. We then get 

*eff = f-a3o-FoGo + a3oGoi^o + — Go + 040-^0 

V 040 / 

+ (022^2 — Osii^oGo + 03iGoi^o + 040-F0 

+ 04oG^ + 04ii^o')e^ (103) 

Go can be eliminated as follows: Vary the 0(e'^) term 
with respect to Go, find the condition 



Ljq — _^o- 



(104) 



If we insert this relation into Eq. (jl03p . the term dis- 
appears after a partial integration and we are left with 



"40 p2 
off - -2-^0 
"30 



2O31O40 — 03004 



1 i?2 



030 



Fo + 040^0 + 022^0 • 

(105) 

Here we have set the formal expansion parameter e = 1 
since it is not needed anymore. The coefficients may be 
simplified by rescaling. 



Foiv) = A/(xy). 



(106) 



2031O40 — 030O41 



MO 



X 



Q30A 

040 



then yields the simpler expression 

^.^=M[{f"f-{.f'f + f + nf] 
with only two residual parameters 

N ~ 040 A^, 

022 1 

— T2- 

040 

Now we focus on the reduced effective action 



(107) 



(108) 



(109) 



(110) 



As we have not been able to solve the Euler-Lagrange 
equation 



(111) 



analytically, we minimize the reduced effective action 
with the Fourier series ansatz 



f{z) =^c„ cos [{2n -f l)qz\ . 



(112) 



Provided we keep only one term in the sum (rimax = 0), 
everything can be worked out analytically with the result 



Co 

q 



1 - 4k 



1 

71' 



(113) 



The (spatially averaged) reduced effective action in this 
approximation is given by 



(V'off) 



1 

96 



(1 - 4.^) 



.^2 



(114) 



The 2nd order phase boundary is obtained from (i/'off) 
and assumes the simple form 



1 

K = —. 

4 



(115) 



The homogeneous, massive solution in the rescaled model 
is characterized by 



g = 0, Co 
and has the reduced action 

V'hom = - 




(116) 



(117) 
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0.6 T 




0.2^ 
0.1- 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 

FIG. 6: Tricritical behavior from GL theory for ^ = 

0.1,0.4,1.0, from left to right. This figure shows how the 
1st order and 2nd order phase boundaries merge in a cusp at 
the tricritical point. 



sufficient for all practical purposes, 

Co = +0.4092971855 

ci = -0.0021185033 

C2 = +0.0000036861 

C3 = -0.0000000064 

q = +0.7064259383 (121) 

Due to the rapid convergence of the Fourier series, the 
lowest order approximation (ri,,„ax = 0) to f{z) is already 
very close to the full result. Likewise, a calculation of the 
spatially averaged effective action, 

(^eff) = -0.0104166667 (n,„ax = 0) 

(V^eff) -0.0104526283 (n,„ax = 3) (122) 

confirms the excellent convergence. 

VII. FULL PHASE DIAGRAM 



The 1st order phase boundary then follows from the con- 
dition (ipcff) = V'hom or 



V6 
4 ■ 



(118) 



Eq. (|118p defines the 1st order phase boundary in the 
tricritical region. Since the final formulae for all coeffi- 
cients and phase boundaries are quite complicated, we 
have collected them in the appendix. These results have 
been used to draw the tricritical behavior for 3 values of 
^ as shown in Fig. [HI 

Truncating the Fourier series (|112p after a single term 
may seem too crude an approximation. Actually, if we 
keep more terms and minimize the effective action numer- 
ically, wc get results which are almost indistinguishable 
on our plot. To illustrate this point, we take the sim- 
pler case where we move along the homogeneous phase 
boundary. According to Eq. (|109p . k ~ Q and the reduced 
effective action (IllOl) becomes 



(119) 



Let us minimize this action using the Fourier ansatz 
pi2p . For n,nax = 0, wc find the analytical result from 
above. 



Co = ^ = 0.4082482904, 
V6 

q = ^ = 0.7071067812. 

V2 



(120) 



For larger values of nmax, the minimization has to be 
done numerically. The following result for rimax = 3 is 



So far, we have discussed only those results about the 
generalized GN model that could be obtained analyti- 
cally, or at least with a minimal numerical effort. For 
the sake of completeness we have also determined the 
full phase diagram with the help of the HF approach for 
a number of values of the parameter ^, interpolating be- 
tween the well-known GN and NJL2 phase diagrams. As 
is clear from the previous sections, for each ^ one needs 
to determine three phase boundaries meeting at the tri- 
critical point: 

• The 2nd order critical line separating massless and 
massive homogeneous phases, identical to the cor- 
responding phase boundary in the original phase 
diagram of the GN model [l^ • This phase bound- 
ary has already been discussed in Sec. IVII and is 
given analytically by Eqs. ()84I85|) . In our case, it 
connects the NJL2 critical point to the critical point 
for a given value of ^, Eq. (|92p . 

• The 2nd order phase boundary separating the soli- 
ton crystal from the massless homogeneous phase 
which can be determined perturbatively (i.e., treat- 
ing the potentials S, P in the Dirac-HF equation 
in 2nd order perturbation theory). The numerical 
work here amounts to one-dimensional numerical 
integrations and solution of transcendental equa- 
tions and can be done easily to any desired accu- 
racy. Moreover, an asymptotic expression for large 
chemical potential will be given in closed analytical 
form. 

• A 1st order phase boundary between crystal phase 
and massive Fermi gas which requires a full numer- 
ical HF calculation. Since the technique has been 
set up previously in a study of the massive NJL2 
model and is described in detail in Ref. Q , we shall 
be very brief here and merely show the final results. 
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Consider the perturbative phase boundary between crys- 
tal and massless Fermi gas first. The calculation is similar 
to the corresponding one for the massive N JL2 model Q , 
except that we may set m = right away. Introducing 
the Fourier components Si , Pi of the HF potentials via 

S{x) = 2Si cos i2pfx) , P{x) = 2Pi sin {2pfx) (123) 

where the Fermi momentum Pf is related to the mean 
fermion density as 



a IT 



(124) 



the single particle energies in 2nd order perturbation the- 
ory read 



En 



(5i-Pi)2 , (Si+Pi^ 



2{p + pj) 2{p-pf) 



ivp > 0) 



= -P - - iw < 0) (125) 

2{P + Pf) 2{p-pf) 

The correction to the single particle contribution of the 
grand canonical potential density is then given by 

/•A/2 

<5*s.p. = P.V. / dp if I + f2 + h + h) (126) 



with 



p{sl + Pi) 


f2 = 


2pfPiSi 
tt{p'^ -pj) 


P{SI + PI) ( 


1 


1 


<P^-p)) V 


-1- qP^p- 


-m) ' 1 + c^(p+m) 


2pfPiSi f 


1 


1 


nip^-p}) [l~ 


h e'3(p- 





(127) 



/4 = 



As in any HF calculation it has to be supplemented by 
the double counting correction, 



-(S'2 + Pf)lnA+^p2 (128) 

TT TT 



Carrying out the principal value integrals involving fi , /2 
analytically, we arrive at the finite expression for the sum 
of dnil) and (fT28ll 

1 f f°° 

6^ = -{Sf + Pi) H2pf) + p_v. / dp(/3 + h) 

(129) 

From here on, we can proceed in the same manner as in 
the previous sections, i.e., we set 

5-^ = MiiSl + 2Mi2SiPi^r M22PI (130) 
and solve the equations 

detA4 = 0, -^detTW = (131) 
dpj 



numerically. Further simplifications occur at large /i 
where the asymptotic behavior of the phase boundary 
can be determined analytically. Once again we take over 
the corresponding formula from the massive NJL2 model 
, merely modifying the double counting correction and 
dropping the S'o(= m) piece. Setting Si = X + y/2, Pi = 
X — y/2, we then get 



2X\ ,. . ,,„,^.2^ .^ , e 



*cff = ^ H^Pf) + H/) - 1) + ^ - f 



— Tdplnfl + e-^V^ 

/^TT Jo V 



(132) 



Minimization with respect to X yields 

X ^ 

41n(4p/) + 2C' 
Minimization with respect to y gives the condition 

1 



(133) 



= 2/ dp 
Jq 

+ \ny 



Vp^TF (1 + c^^^^^) 

ein(4p/) 



e + 21n(4py-)' 
Expanding the integral in (|134p for small y [1 
Cln(4j3/) 



(134) 



= Iny 



C + 21n(4p/) 



C + 0(y^), (135) 



the asymptotic form of the phase boundary is finally 
given by the expression (/^ « py). 



-K 



K 



ein(4/i) 
^ + 2 ln(4/i) ■ 



(136) 



X in Eq. (|133l) interpolates between (NJL2) and y/2 
(GN) for ^ = O...00. Likewise, Tcrit smoothly interpolates 
between the known results for the NJL2 and GN model, 
respectively. 

In Fig. [7] we show by way of example the perturba- 
tive phase boundary at ^ = 1.2, together with the NJL2 
(^ = 0) and GN (^ -> 00) model phase boundaries. The 
asymptotic expression p36p is shown as the dashed curve 
and only deviates from the full result below /t « 1. Fig. [5] 
represents a 3d plot of the perturbative phase boundary 
for 10 values of ^ ranging from to 10. The thick line 
is the tricritcial curve. We have also drawn asymptotic 
behavior according to Eq. (|136p for 3 moderate values of 
ji to demonstrate how well this simple formula catches 
the perturbative critical sheet for all values of ^, starting 
from /i « 1. 

Still missing in Fig.|S]is the critical sheet separating the 
crystal from the massive Fermi gas. We recall that this 
phase transition is of 2nd order in the GN model, non- 
existing in the massless NJL2 model and of 1st order 
in the massive NJL2 model. Wc find that it is of 1st 
order in the generalized GN model for all values of ^, so 
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GN 



0.5 1 1.5 2 

FIG. 7: Perturbative 2nd order phase boundary separating 
the crystal from the chiraUy restored homogeneous phase at 
^ — 1.2. Also shown are the corresponding critical lines for 
the NJL2 model = 0) and the GN model = 00). Dashed 
curve: asymptotic expression, Eq. (|136|l . The open circles are 
the tricritical points for all 3 cases. 




FIG. 8: Like Fig.El 3d plot for several values of ^ = 0, 0.1, 
0.2, 0.4, 0.8, 1.2, 2.0, 3.0, 5.0, 10.0). Fat curve: Tricritical line. 
Three curves at constant /i (/i = 1.0, 1.5, 2.0): asymptotic 
expression, Eq. (|136p . 



that apparently the phase transition becomes continuous 
only in the GN limit ^ — > 00. Hence there is no way 
of determining the critical sheet perturbatively and we 
need a full thermal HF calculation. Fortunately, this can 
be done using the techniques which have recently been 
developed for the massive NJL2 model [1]. As a matter 
of fact, all what is needed is a trivial modification of the 
double counting correction. We therefore refer to Ref. Q 
for more technical details and immediately pass on to the 
results. 



Let us first consider the 1st order critical line at T = 0, 
i.e., the baseline of the 1st order critical sheet in a 3d plot. 
This is closely related to the baryon mass discussed in 
Sec. IIVI near the chiral limit. Since we are not restricted 
to small ^ values in the numerical HF calculation, we 
can now get complementary information to the one of 
Sec. IIVI and complete the picture about baryons in the 
generalized GN model. Fig. IH] shows the phase bound- 
ary at zero temperature in the (^, /z) plane (the actual 
calculation was done at T = 0.05, but this makes no 
difference). Since baryon number is 1/2 in our model, 
the critical chemical potential has to be identified with 
twice the baryon mass (divided by N) here. The reason 
is the following: The critical chemical potential at T = 
is the amount of energy needed to add a fermion to the 
vacuum. If the kink-like baryon has mass Mb and car- 
ries N/2 fermions, we get Hcrit = 2Mb/N. The curve in 
Fig. ini interpolates between the massless baryons of the 
NJL2 model and twice the mass of the kink in the GN 
model, Mb/N = I/tt. As shown in Fig. [TUl at small val- 
ues of ^ the numerical HF results match nicely onto the 
derivative expansion, a welcome test of both the analyt- 
ical and numerical approaches. From the HF calculation 
at the phase boundary we can also extract the shape 
of the self-consistent potentials for a single baryon, now 
for arbitrary values of ^. A typical example is shown in 
Fig. [m for the case ^ = 2. The scalar potential has kink 
shape at all ^, going over into the GN model kink in the 
limit ^ — > 00. The pseudoscalar potential is bell shaped 
and gets more and more suppressed with increasing ^. 
This is of course just the effect of the double counting 
correction term (|33)) where ^ acts like a Lagrange multi- 
plier for P, quenching it completely in the limit ^ ^ 00. 
The other limit, — > 0, has already been discussed before 
in Sec. IIVI in terms of the sine- Gordon kink with scalar 
and pseudoscalar potentials of the same amplitude. 

Finally, we come to the full phase diagram as a func- 
tion of ^, /i, T, including the numerically determined 1st 
order sheet. It is shown in Fig. [12] and Fig. [13] under 2 
different viewing angles for the sake of clarity. As ex- 
plained in more detail in Ref. 0, the phase boundary 
is determined by performing the HF calculation along 
a trajectory crossing the critical line and comparing the 
grand canonical potential of the massive Fermi gas to the 
one of the soliton crystal. As we know the exact location 
of the tricritical point in the present case, we are even 
in a somewhat better position here than in the previous 
study of the massive NJL2 model. 



VIII. SUMMARY AND CONCLUSIONS 

In this paper, we have studied a generalization of the 
GN model with two different (scalar and pseudoscalar) 
coupling constants. This equips us with an "interpolat- 
ing field theory" between the well-studied massless GN 
and NJL2 models in a way which always keeps the dis- 
crete Z2 chiral symmetry intact. The continuous chiral 
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FIG. 9: First order phase boundary separating the crystal 
from the massive Fermi gas phases at T = in the generalized 
GN model. The vertical axis may be interpreted either as 
critical chemical potential or twice the baryon mass, due to 
fractional baryon number 1/2 in this model. The straight 
line shows the asymptotic value 2/n taken from the standard 
GN model. Numerical calculations performed for a few extra 
points (C = 0.3, 0.6, 1.0, 1.6, 2.5, 4.0, 6.5, 8.0) in addition to the 
values mentioned in the caption of Fig. [S] 



0.7 1 



0.6 




" 0.2 0.4 0.6 0.8 1 

FIG. 10: Like Fig. |9] but blowing up the region of small ^ 
to check the consistency between the derivative expansion of 
Sec. IIVI (lower curve: NNNLO, upper curve: NNLO) and the 
numerical HF calculation (circles). 



symmetry of the NJL2 model is only recovered for equal 
coupling constants, so that we now break chiral symme- 
try (explicitly) in a quite different manner than via the 
usual fermion mass term. Our motivation was primar- 
ily to get further insights into the solitonic aspects of 
4-fermion theories in 1-|-1 dimensions which have been 
investigated intensely in recent years. 

The first insight is the emergence of the dimension- 




i''''2''''3''''4''''5''''6''''7 
X 



FIG. 11: Example of numerical baryon HF potentials at ^ = 
2.0. For larger values of ^, P decreases and S approaches the 
GN kink (not shown). 




FIG. 12: Like Fig. [HI but including 1st order phase boundaries 
separating the crystal from the massive Fermi gas. The curve 
drawn at ^ = 2 is the asymptotic expectation according to 
Eq. (|136p , the base line at T = is taken from Fig. (9] 



less parameter ^ during the process of regularization and 
renormalization, in addition to the familiar fermion mass. 
The basic relations, Eqs. (|16p . which generalize the stan- 
dard gap equation remove all divergences encountered in 
subsequent applications, both in the treatment of bound 
states (mesons, baryons) and in the thermodynamics of 
the model. The parameter ^ plays a role analogous to 
the "confinement parameter" 7 in massive GN models. 
This is particularly striking in the RPA approach to the 
pseudoscalar fermion-antifermion bound and scattering 
states, where the results for the massive NJL2 model and 
the generalized GN model become identical if we replace 
7 by ^. The qualitative effect of ^ on the HF calcula- 
tions at zero and finite temperature is very easy to un- 
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FIG. 13: Like Fig. 1121 but different orientation for better 
visibility 



dcrstand. It only enters in the double counting correction 
to energy or thermodynamic potential as an extra term 
~ ^ / dxP^. Hence it acts like a Lagrange multiplier for 
the pseudoscalar potential, leading to a complete quench- 
ing of P in the GN limit ^ — >■ oo. Thus ^ may be thought 
of as a "chiral quenching parameter" responsible for the 
transition from complex condensates living on the chiral 
circle in the NJL2 model to the purely real condensates 
of the GN model. 

As far as baryon structure is concerned, the most in- 
teresting result is perhaps the fact that the new baryons 
interpolate between the kink of the GN model and the 
massless baryon of the NJL2 model, always carrying frac- 
tional baryon number 1/2. This is certainly a conse- 
quence of the fact that the generalized GN model still 
has a discrete chiral symmetry. Indeed in the massive 
NJL2 model, chiral symmetry is explicitly broken by the 
mass term without a residual Z2 symmetry and one finds 
baryons with integer baryon number 1. This new kind 
of chiral kink is different from all known multi-fermion 
bound states in the GN model family and has been deter- 
mined analytically for small ^ and numerically for large 

The phase diagrams of the NJL2 and GN model look 
very different, so that we were curious to see how our 
theory would manage to interpolate between these two 
pictures. This can now be answered most clearly by 
the study of the tricritical behavior near the chiral limit, 
largely analytically owing to the GL approach. The rel- 
evant picture is Fig. [SJ showing a kind of "morphing" 
from GN-type behavior to the NJL2 phase diagram with 
its single straight line phase boundary. Together with 
the numerical HF calculation, we are now confident that 
the solitonic crystal phase is separated from the mass- 
less (massive) Fermi gas by a 2nd (1st) order transition, 
respectively. This was not clear a priori, since the transi- 
tion from the crystal to the massive homogeneous phase 



is continuous in the GN model and doesn't even exist in 
the NJL2 model. Our interpolated phase diagram also 
looks qualitatively different from the one of the massive 
NJL2 model which has only 2 phases (no massless phase 
due to explicit breaking of the Z2 symmetry), and where 
the opening angle between the 2 phase boundaries at the 
tricritical point was tt rather than 0. 

Initially, we had hoped that the generalized GN model 
can be solved analytically for arbitrary ^, since this is 
what happens at the "endpoints" ^ = (NJL2) and 
^ = 00 (GN). However, this does not seem to be the 
case. In this situation, the fact that our toolbox also 
contains the numerical HF method has turned out to be 
a definite advantage. A combination of analytical calcu- 
lations and a numerical approach gives us confidence that 
we have solved and understood the model in the large N 
limit fairly well. The most serious limitation at present 
is the fact that our techniques are tailored to point-like 
4-fermion interactions and cannot deal with gauge theo- 
ries in a systematic fashion. This is unfortunate in view 
of the interesting features of, e.g., the 't Hooft model [2l| 
where more analytical insights into the early (llj and very 
recent [13, H^l numerical HF calculations on the lattice 
would be welcome. 
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Appendix: Details of the Ginzburg-Landau 
approach of Sec. IVII 



Here wc collect the detailed formulae used in preparing 
Fig.|6]in Sec. IVII We first list the coefficients of the Taylor 
expansion (|98p . Using the notation 



*(")(zO=^n + i/n, zt = T: + 7^, (137) 

z Zn 
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one finds 
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The scale parameters A, x from Eq. (|107p and the residual 
parameter k in the effective action piOp then become, 



A = 



8^(2i?2 + /1/3) 



nRi 



X 



6ATT^Tth{2Rl + hh) 



VlRl 



Rl{na\R2-ll)-vm^) 
ATiaH2Rl + hhf 



(140) 



2nd order phase boundary in local coordinates a, rj, see 

Eq. (nnni), 



V = 



{R2-II)- 



{2Rl + hh 

m 



a^. (141) 



1st order phase boundary, 



.-^((i^.-/?) + (2 + ^/6)M±ilM!V^ 



Rl 



(142) 

These critical lines can easily be rotated back to the orig- 
inal coordinates, see Fig. [6] for some results. 
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